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Abstract 



^ ! We study possible so far unaccounted isospin breaking effects in the relation between the 

pion form factor as determined in e^e~ experiments and the corresponding quantity ob- 
■ tained after accounting for known isospin breaking effects by an isospin rotation from 

^ ■ the T-decay spectra. In fact the observed 10% discrepancy in the respective pion form 

factors may be explained by the isospin breaking which is due to the difference between 
masses and widths of the charged and neutral p mesons. Since the hadronic contribu- 
tion to the muon anomalous magnetic moment can be calculated directly in terms of 
p i' the e'^e^ -data the corresponding evaluation seems to be more reliable. Our estimate is 

p • Qhad(i) _ (^gg4_g -|- 3_g^ X 10~^°. The r-data are useful at the presently aimed level of 

D . accuracy only after appropriate input from theory. 

: 1 Introduction 

S ■ 

The most precise measurement of the low energy pion form factor in e~^e -annihilation 
experiments is from the CMD-2 collaboration. The updated results for the process e~^e^ 
p vr+vr" have just been published fl]. The update appeared necessary due to an 
overestimate of the integrated luminosity in previous analyses. The latter was published 
in 2002 0- A more progressive error estimate (improving on radiative corrections, in 
particular) allowed a reduction of the systematic error from 1.4% to 0.6 % . 

Since 1997 precise r-spectral functions became available |S1 IH E] which, to the extent 
that flavor SU{2)f in the light hadron sector is a symmetry, allows to obtain the iso- 
vector part of the e"'"e~-cross section (Oj. In this way r data may help to substantially 
improve our knowledge of [^^^^(s), which is an important input for the evaluations of 
the hadronic vacuum polarization contributions to the anomalous magnetic moment of 
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the muon and of the effective fine structure constant aem(j^z) an important input for 
LEP/SLC precision physics (see e.g. [7]). The idea to use the r spectral data to improve 
the evaluation of the hadronic contributions a^'^ and Aa^^'^ was pioneered in jH]. 

With increasing precision of the low energy data it more and more turned out that 
we are confronted with a serious obstacle to further progress: in the region just above 
the a;-resonance, the isospin rotated r-data, corrected for the known isospin violating 
effects jS], do not agree with the e+e~-data at the 10% level jTU]. Before the origin of this 
discrepancy is found it will be hard to make progress in pinning further down theoretical 
uncertainties in the predictions for and aem(^z)- 

In this context isospin breaking effects in the relationship between the r- and the 
e"'"e~-data have been extensively investigated in jH]. One point which in our opinion has 
not been satisfactorily clarified is the role of the isospin breaking effects in the charged 
vs. neutral p line-shape, which must manifest themselves in mp± — rripO and Fpi — TpO. 
Looking at the particle data tables [TT], there is no established non-zero result as yet. 
Earlier statements about the problem in |HJ 13 CHI adopted essentially the PDG estimate 
mp± —rripO = Oil MeV. There are theoretical arguments about why this mass difference is 
expected to be very small: the usual argument assumes Am^ ^ Am^ via a sum rule, which 

then yields nip- — uipo ~ ^^UhL ^ 0.014 MeV ! This is based, however, on an assumption 

which need not be true for the mass definition adopted in recent p-line shape analyses^. 
If fact a more recent analysis of the CMD-2 (before the last update) and the ALEPH and 
CLEO data yielded mp± - rupo = 2.6 ± 0.8 MeV and Tp± - TpO = 3.1 ± 1.7 MeV 
where the uncertainty is our estimate. The corresponding isospin corrections may still 
look too small to account fully for the observed discrepancy in the spectral functions but 
they clearly point towards a substantial reduction of the problem. Our strategy therefore 
here is a different one. Our hypothesis is that as a leading effect the discrepancy very 
likely is due to the isospin breaking by the charged vs. neutral p-meson parameters. A 
similar but subleading contribution is expected to come from possible isospin violations 
in the respective backgrounds (encoded usually by the p',p" contributions). Since the 
fit formulae adopted, like the Gounaris-Sakurai formula UHl; are far from being based 
on first principles we should not trust to much in the fitting procedures based on them. 
E.g., usually just a set of resonances p, u, p', p" is included but we have no idea about the 
background (continuum) which also should be included somehow. We also would like to 
advocate that, as the level of accuracy of the present discussion advances, in future one 
should compile charged and neutral p data separately. 

Whether the observed discrepancy is an experimental problem, or just a so far un- 
derestimated isospin breaking effect will also be settled, hopefully, by new results for 
hadronic e~^e~ cross-sections which are under way from KLOE, BABAR and BELLE. 
These experiments, running at fixed energies, are able to perform measurements via the 
radiative return method jTTl UHl HH] • Results presented recently by KLOE seem to agree 

^For a more detailed theoretical estimate see jj^l and references therein (see also ^| ^J). Mass 
and width of an unstable particle, and in particular of the p, depend on the precise definition. In that 
sense they are pseudo-observables which depend on theoretical input. What we need here is a consistent 
prescription to extract them from the experimental data (p-line shape) and make sure that we compare 
comparable quantities (with respect to the handling of vacuum polarization effects, final state radiation, 
energy dependence of width, background etc.). 
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very well with the final CMD-2 e^e -data. 



2 The T vs. e+e problem 

The iso-vector part of a{e^e^ — >■ hadrons) may be calculated by an isospin rotation from 
r-decay spectra, to the extent that the so-called conserved vector current (CVC) would 
be really conserved (which it is not, see below). The relation may be derived by comparing 
the relevant lowest order diagrams which for the e^e~ case translates into 



(0) - t + - + -\ ^vra^ 



and for the r case into 
-— r TT ttV^ = — — 1 -] 1 + — v_{s) 2 

where \Vud\ = 0.9752 ± 0.0007 ^T] denotes the CKM weak mixing matrix element and 
'S'EW(new) = 1-0233 ± 0.0006 [S'Ew(oid) = 1.0194] accounts for electroweak radiative correc- 
tions 1^ 1221 12ni ini ^] • The spectral functions are obtained from the corresponding 
invariant mass distributions. The B{i)^s are branching ratios. SU(2) symmetry (CVC) 
would imply 

v.{s) = vo{s) . (3) 
The spectral functions Vi{s) are related to the pion form factors F^{s) by 

v^is) = ^\FM' ; (^ = 0,-) (4) 

izvr 

where Pf{s) is the pion velocity. The difference in phase space of the pion pairs gives rise 
to the relative factor (3^.^0/(31-^+ UM- 

It is important to check what precisely the experimental data in each case represent. 
In CMD-2 e'^e' measurements we exclude final state radiation (FSR) as well as vacuum 
polarization effects^ which are not included in the r data in first place. FSR as far as 
included by the measurement (virtual and soft real radiation) has been subtracted together 
with the initial state photon radiation. Hard FSR photons were rejected to a large extent. 
In our analysis Ft^{s) obtained from CMD-2 is the undressed (from VP and FSR) pion 
form-factor (see e.g. |2H])- The available r decay spectra all include photon radiation (no 
subtractions were made), which hence has to be subtracted a posteriori (the correction 
factor Gem below), while photon vacuum polarization effects play no role (i.e., are not 
included). This is because the r decay as a charged current (CC) process proceeds by 
the heavy W exchange, which makes it an effective four-fermion interaction with Fermi 
constant Gp as a coupling in place of a{s). In contrast to a the Fermi coupling Gp is 
not running up to LEP energy scales. Electroweak short distance corrections (hadronic 

^Tab. 1 of fj lists IFttP, which includes VP but not FSR, as well as cr^^^^-, which includes FSR but 
not VP. VP and FSR are separately known quantities which we may add and subtract according to our 
needs. 
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relative to leptonic channel) give rise to the correction factor Sew = 1 + (^ew? which is 
dominated by a leading large logarithm (l + (a/7r) \ia{Mz/mr)) which should be resummed 
using the renormalization group [20] • Note that the overall coupling drops out from the 
ratios in This also makes it evident that the subtraction of the large and strongly 
energy dependent vacuum polarization effects (see e.g. Fig. 1 in [20]) necessary for the 
e"'"e~-data, which seems to worsen the e^e~ vs. r problem, was properly treated in 
previous analyses. 
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Figure 1: Isospin corrections applied to the r data: left the corrections [9., with tib 
defined in ([2]), and right the barely visible effect on the ALEPH data. 



Before a precise comparison via ([HJ) is possible all kind of isospin breaking effects have 
to be taken into account. As mentioned earlier, this has been investigated carefully in jOj 
for the most relevant tttt channel. Accordingly, we may write the corrected version of ^ 
(see |9j for details) in the form 



(0) 



with 



Kr 



Gl 



a. 



\Vud?ml 



Kris) 



r[7] 



ds 



X 



Rmjs) 
Sew 



3847r3 

and the isospin breaking correction 



1 - 



1 + 2 



Rmis) 



1 (3^. 
Gem{s) pi 



Fyis) 



~3s 



(5) 



(6) 



(7) 



The factor G'em(s), displayed in Fig. [H includes the QED corrections to — > z/t-tt^tt^ 
decay with virtual plus real soft and hard (integrated over all phase space) photon radi- 
ation calculated in scalar QED, except from the short distance term, which is calculated 
for the corresponding quark production process and included conventionally in Sew as 
the leading logarithm, mentioned before. Originating from (H)), /?^_.^+//3^_^o (see Fig. [T]) 
is a phase space correction due to the tt^ — 7r° mass difference. Fv{s) = F^{s) is the 
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NC vector current form factor, which exhibits besides the 1 = 1 part an / = con- 
tribution. The latter p — u mixing term is due to the SU(2) breaking (m^ — m„ mass 
difference). Finally, f+{s) = F~ is the CC 1 = 1 vector form factor. One of the leading 
isospin breaking effects is the p — u mixing correction included in The transition 

— > \Fy~^\s)\'^ ~ |/+(s)p is illustrated in Fig. E| The form-factor corrections, 
in principle, also should include the electromagnetic shifts in the masses and the widths 
of the p's. Up to this last mentioned effect, which was considered to be small, all the 
corrections were applied in |IOj but were not able to eliminate the observed discrepancy 
between and Vo{s) (see [TU] for details and Fig. 01 below). 

In fact, possible isospin breaking corrections due to different electromagnetic shifts of 
masses and widths of the neutral and charged p-mesons (remember that for other quark 
bound states the vr's we have m^± — 771^^0 = 4.5935 ± 0.0005 MeV), respectively, have 
been mentioned or were very briefly discussed only in |S1 IHl and they might have 
been underestimated so far. Such isospin violating mass and width differences are not 
established experimentally, not even the sign. We therefore ask the question whether 
applying an isospin breaking correction which accounts for that could resolve the puzzle 
of the above mentioned discrepancy. 

What we do is the following: we take the CMD-2 data and subtract off the VP and 
the cu-contribution. The latter / = part enters via p — uj mixing, which is a consequence 
of isospin violation due to the mass difference rriu — md of the light quarks. To this end 
we may take the Gounaris-Sakurai kind parameterization (see [THl |2] for details) (FSR 
not included; BW=Breit-Wigner; BW'~^^=Gounaris-Sakurai modified Breit-Wigner ) 



FJs) 



BW;rf,,o)(^) ■ (1 + 5i7^BW^(.) + P BWf,,,,^{s) + 7 BW%,,^{s) 



l + /? + 7 

of the CMD-2 data and set the mixing parameter 6 = 0. In this way we obtain the 
iso-vector part of the square of the pion form factor |F^|^g+'^3_-)(s) displayed in Fig. |21 To 
the r version of the pion form factor, following from Q and we perform the isospin 
breaking corrections 

f3n-TT+ 'S'EW(old) 



riBis) = ^ (9) 



5 



1.30^ 

1.25 
s 1.20 
5 1. 15 
- 1. 10 

1.05 

LJ 

-,1.00 

^ 0.95 
0.90 




ratio of the fits 

Rel. [2/1 




1000 1100 



600 650 700 750 800 850 
E (MeV) 



950 1000 



Figure 3: Left: The ratio between r-data sets from ALEPH, OPAL and CLEO and the 
1 = 1 part of the CMD-2 fit of the e"'"e~-data. The curves which should guide the eye are 
fits of the ratios using 8th order Tschebycheff polynomials. Right: Ratio of the ALEPH 
vs. CMD-2 fits differing by mass and width only (see Tab. Q). By the isospin violation 
correction {mp~,Tp-) (mpOjPpo) of the r-data this ratio becomes trivially equal to 
unity. 



with Gem{s) from 0. The factor 5'EW(oid) /5'EW(ncw) corrects for some previously missing 
corrections CHI- The such obtained corrected'^ pion form factor |-F7r|^/)^^('5) is to 
be compared with \Ft^\'^J-^J-^{s) . The ratio shows the unexpected large deviations from 
unity (see Fig. El). While the ALEPH and CLEO data clearly exhibit the structure as 
expected from an increase of mass and width of the p the OPAL data show a different 
form of the spectrum. The problem with the OPAL data originates from the fact that 
in the neighborhood of the p peak the cross section apparently is too low. Since the 
distribution is normalized to the very precisely known total branching fraction, the tails 
of the resonance get enhanced, which leads to the structure actually seen (the apparent 
width gets enhanced). Of course the data points of the spectrum have imposed strong 
error correlations via the normalization to the integral rate. In the figure only the diagonal 
elements of the covariance matrix are visualized. 

For a comparison of the e~^e~ with the r data we have to perform fits in a common 
energy range. It is important to fix the "p-background" appropriately. The latter is 
represented in the GS-parametrization by the higher resonances p and p and in order to 
fit the corresponding parameters we have to extend the fit range towards higher energies. 
To this end we are including beyond the CMD-2 limit the available vr+vr" data from other 
experiments up to 1.5 GeV above which the r data get of low quality. We finally decided 
to adopt the PDG values for mp> = 1465 ±25, Fp/ = 310 ±60 and mp" = 1700 ±20, Fp// = 
240 ± 60 for the masses and widths and fit the complex admixing coefficients f3 and 7. 
For \f3\ = 0.12 ± 0.005, (pf^ = 160 ± 5 and I7I = 0.023 ± 0.005, cp^ = we get reasonable 
fits for all data sets. In particular the widths depend substantially on the precise values 
of these parameters. Fortunately, the parameter shifts of interest turn out to be rather 
stable. In other words, while there are different possible parametrizations of comparable 
quality for the individual data sets the systematic shift between e~^e~ and r data is there 

•^The velocity factor correction of course only applies when, as frequently has been done, the wrong 
velocity was used in Q in the extraction of the charged channel form factor. 
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T ALEPH 


T CLEO 


r OPAL 


e+e- CMD-2 


rUpO 
Tpo 
rup- 

S 


775.52 ± 0.49l[J:i 
149.40 ± 0.68l[J:J[| 
1.39 


775.01 ± 0.36t[j:i^ 
149.00 ± 0A9tlil 
1.35 


777.34 ± 1.2ll[;:i 
153.91 ± l.62l:\fo 
0.62 


772.95 ± 0.56 ± 0.12 

147.93 ± o.rollil 

1.28 



Table 1: Results of fits to the isospin breaking corrected pion form factors squared for 
r (ALEPH, CLEO and OPAL) and e+e^ (CMD-2) data. The Gounaris-Sakurai parame- 
terization of the p line shape is utilized. Masses and widths in MeV. S = Jx^ / {f^ ~~ 1) 



independently of how we parametrize the data. 

We may fit now \Fj^\'^^^^^{s) with the Gounaris-Sakurai formula (jHI) with no uj term, 
i.e., with 5 = 0, in order to obtain mp± and Tp±. We would like to emphasize that it is 
important to "zoom-in the p" appropriately in determining m^i — iripo and T p± — PpO, 
i.e., we have to perform the fits at fixed background (i.e., besides and Tp all other 
parameters in the GS-formula are held fixed) in the p dominated region between 610.5 
MeV and 961.5 MeV (CMD-2 range). 

This simple leading effect analysis yields the results in Tab. As can be seen the 
r data give consistently larger values for both mass and width of the charged p. The 
evidence is far from impressive, between ALEPH and CMD-2 we have Am^ = 2.6 ± 0.8 
MeV and AFp = 1.5 ± 1.0 MeV. The two-parameter fits are not of good quality and 
a more elaborate analysis would be needed to come to more precise conclusions. An 
analysis based on more recent preliminary ALEPH data, yields the slightly larger values 
Amp = 3.1 ± 0.9 MeV and APp = 2.3 ± 1.6 MeV |2Z|; which are consistent however with 
our findings. 

Note that our values of nipO and PpO differ substantially from the values 775.65 ± 0.64 ± 
0.5 and 143.85± 1.33±0.80, respectively, given by CMD-2. The latter have been obtained 
by fitting the physical vrvr cross-section (before subtracting the VP) and including the uj 
and the p' with mixing parameters \5\ = (1.57 ± 0.16) x 10~^, ips = 13.3° ± 3.7° and 
|/?| = 0.0695 ± 0.0053, (pi3 = 180°. A p" was not included. The parameters determined by 
CMD-2, cannot be directly compared with corresponding parameters obtained by fits of 
the T data, because the latter do not include VP effects. In order to have a common basis 
we subtract the VP from the e'^e~ data (of course, alternatively, we could supplement 
the r data with the VP [by replacing a in K^j (0) by a{s)] in order to obtain dressed 
(physical) parameters which are the ones usually listed in the particle data tables). The 
subtraction of the VP lowers the mass by about 1 MeV and increases the width by about 
1.3 MeV. The additional changes are due to the inclusion of the p" and the changes in the 
other "background parameters" . If we would utilize the CMD-2 background parameters 
as a common background parametrization we would not be able to get acceptably good 
fits for the r data sets. 

''For the e^e" channel we fit the data after subtraction of the 1 = part. Including the uj in the fit 
yields practically the same result. 
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Now we assume that the systematic deviations seen in the p parameters include 
electromagnetic isospin breaking which we have to correct for. We now may ask two 
questions. The first is: how does the test-ratio of Fig. |21 look like if we replace mp± and 
T p± in the r-data fit by the more appropriate rupo and T po ? The second is: what mass 
and width do we get if we fit them in the r-data parameterization such that the test-ratio 
comes out to be unity within errors? Not too surprisingly we find them close to the ones 
given in Tab. Qfor CMD-2: It makes the central value of the ratio unity within 0.1 % ! 
Uncertainties may be obtained from the ones in the parametrizations. Of course keeping 
the background fixed the result looks pretty trivial. In fact fitting all parameters of the 
GS formula simultaneously in the much wider range of r data, as has been performed 
in ^21; yields results which look very similar to ours. The parameters obtained are very 
strongly correlated and all of them may be affected by some isospin breaking effects. A 
much more elaborate analysis would be necessary to actually establish tight experimental 
values for possible isospin breakings in these parameters. 

We conclude that the r vs. e^e~ discrepancy very likely is an isospin breaking effect 
which has not been accounted for correctly in previous analyses. This also would establish 
a significant difference for mp± — rripo and Fpi — TpO. Of course, what it means is that the 
r data cannot be utilized to calculate a^^'^ without reference to the e+e~ data. Also, since 
now substantially correlated, the inclusion of the r-data is much less straightforward. 
The question is how much they still can contribute to reduce the uncertainties in the 
evaluation of a^^'^. This also makes it very likely that the e"''e~-data based evaluations 
are the more trustworthy ones. After the correction in the normalization of the CMD-2 
data we get the leading hadronic contribution to the anomalous magnetic moment of the 
muon. We now obtain 

aJ^^W = (694.8 ± 8.6) X 10"^° [e+e" - data based]. (10) 

With this estimate we get 

af'' = (11 659 179.4 ± 8.6had ± 3.5lbl ± 0.4qed+ew) x 10"^° (11) 

which compares to the most recent experimental result [2H] 

a^fP = (11 659 203 ± 8) x 10^^° . (12) 

The "discrepancy" \a^f^'^ — a'^p'^l = (23.6±12.3) x 10^^° corresponds to a deviation of about 
1.9 cr. For other recent estimates we refer to pUj I29j. 

3 Summary and Conclusion 

Since recently we have in each case two reasonably consistent sets of data: the ALEPH 
and CLEO r-data sets on the one hand and the CMD-2 and KLOE (still preliminary) 
e"'"e~-data sets on the other hand. The r-data samples are about 10% higher than the 
e~^e~ ones in the tail above the p. This can be clearly seen in Fig. |21 Assuming that 
the experiments are essentially correct we think that a 10% increase in a resonance tail 
can easily be attributed to a 0.5% increase in the energy scale. Since it is very unlikely a 
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problem of energy calibration, the only explanation remains that the resonance parameters 
must be different in the charged and the neutral channel. Our estimated shifts in the p 
parameters account for about 6.8% (or 8.1% with the estimates given in [TTj) in the cross 
section (see Fig. El). The remaining "discrepancy" is likely due to corresponding shifts in 
the other GS fit-parameters. 

Our analysis shows that to a large extent we may understand the e^e~ vs. r discrep- 
ancy as an isospin breaking effect coming from the fact that mass and width of charged 
and neutral p-mesons, as naively expected, are different and thus that the r-data must 
be mapped to the neutral channel parameters before they can be utilized for the eval- 
uation of a^'^ in addition to the e"'"e~-data. We thus assume that a main part of the 
problem is due to additional isospin breaking effects and not primarily an experimental 
one. Of course there are also experimental difficulties which hopefully will be resolved by 
forthcoming experiments. 

The fact that the p-mass difference is found to be of size comparable (maybe half of it) 
to the well-established pion mass difference m^ri — m^o ~ 4.6 MeV seems not so unlikely 
because the corresponding bound states, apart from the spin orientation, have the same 
quark content (p^ vs. vr^ on the one hand and p° vs. 7r° on the other hand)^. In the 
charged case the electric potential is repulsive while in the neutral case it is attractive 
and contributes to lower the mass, irrespective of the spin orientation of the constituents. 
We do not think that the Goldstone-boson nature of the pions, which derives from the 
properties of the strong interactions, the spontaneous breaking of the chiral symmetry, 
necessarily makes such a manifest isospin breaking effects completely different for the 
vector particles and the scalars. 

Since we do not have an independent evaluation of the charged and neutral p-meson 
parameters, the isospin correction needed in order for the r-data to be useful for the 
evaluation of the hadronic contribution to the muon anomalous magnetic moment cannot 
be performed at sufficient precision at the moment. Nevertheless, the r-data still provide 
important cross checks and last but not least Ref. jS] triggered a discussion which forced 
all parties to check more carefully what they have done. One impact was that also the 
results in the e"'"e~-channel had to be corrected. 

The main point is that there must be such effects which were not correctly treated 
so far. Maybe the available data do not allow us to pin down a solid value for the mass 
and width differences (all GS parameters in fact must be subject to isospin breaking and 
the data may not suffice to come to a definite conclusion). We think the main point is 
that the derivatives of with respect to rup and Vp are rather large and hence no stable 
result can be given if the isospin breaking in these parameters is not known with sufficient 
precision. 

Note that in spite of the fact that the dominating p-peak is shifted downwards, due 
to the correction which we have to apply to the r-data, the weighted a^-integral does 
not increase. It rather decreases, because the width also substantially decreases by the 
correction and actually over-compensates the effect of the shift in the mass. 

^Remember that, for example, para- (spin s = 0) and ortho- (s = 1) positronium have almost the same 
binding energy. For the pion and the p the average distance of the quarks is determined by the strong 
interactions, and in principle could be different. However, the charge radius of the p is of similar size as 
the one of the pion and thus also the electromagnetic effects are expected to be of similar magnitude. 
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Our conclusion: very likely we are back to one prediction for which is the e~^e~- 
based value at about 2a below the experimental result! Unfortunately, at present, we do 
not have a precise enough understanding of the isospin violations to be able to utilize the 
r-data for the evaluation of the hadronic contribution to (7 — 2 of the muon. 
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